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Abstract 

We develop a sufficient condition for the least-squares measurement (LSM), or the square- 
root measurement, to minimize the probability of a detection error when distinguishing between 
a collection of mixed quantum states. Using this condition we derive the optimal measurement 
for state sets with a broad class of symmetries. 

We first consider geometrically uniform (GU) state sets with a possibly nonabelian generating 
group, and show that if the generator satisfies a certain constraint, then the LSM is optimal. 
In particular, for pure-state GU ensembles the LSM is shown to be optimal. For arbitrary GU 
state sets we show that the optimal measurement operators are GU with generator that can be 
computed very efficiently in polynomial timc, within any desired accuracy. 

We then consider compound GU (CGU) state sets which consist of subsets that are GU. 
When the generators satisfy a certain constraint, the LSM is again optimal. For arbitrary CGU 
state sets the optimal measurement operators are shown to be CGU with generators that can 
be computed efficiently in polynomial time. 
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1 Introduction 



In a quantum detection problem a transmitter conveys classical information to a receiver using a 
quantum-mechanical channel. Each message corresponds to a preparation of the quantum channel 
in an associated quantum state represented by a density operator, drawn from a collection of known 
states. To detect the information, the receiver subjects the channel to a quantum measurement. 
Our problem is to construct a measurement that minimizes the probability of a detection error. 

We consider a quantum state ensemble consisting of m density operators {pi, 1 < i < m} on an 
n-dimensional complex Hilbert space Tí, with prior probabilities {pi > 0, 1 < i < m}. A density 
operator p is a positive semidefinite (PSD) Hermitian operator with Tr(p) = 1; we write p > to 
indicate p is PSD. A mixed state ensemble is one in which at least one of the density operators 
Pi has rank larger than one. A pure-state ensemble is one in which each density operator pi is a 
rank-one projector \<pi){(pi\, where the vectors \(f>i), though evidently normalized to unit length, are 
not necessarily orthogonal. 

For our measurement we consider general positive operator- valued measures Q ^]. Necessary 
and sufhcient conditions for an optimum measurement minimizing the probability of a detection 
error have been derived (3|, ^ ||. However, in general, obtaining a closed-form analytical expression 
for the optimal measurement directly from these conditions is a dimcult and unsolved problem. 
Iterative algorithms minimizing the probability of a detection error have been proposed in |6[ [5| . 

There are some particular cases in which the solution to the quantum detection problem is 
known explicitly |ï|, [?|, ||, ||, 10 1. Ban et al. Sj derive the solution for a pure-state ensemble 



consisting of density operators pi = \4>í){4>í\ where the vectors \4>i) form a cyclic set, í.e., the 
vectors are generated by a cyclic group of unitary matrices using a single generating vector. The 
optimal measurement coincides with the least-squares measurement (LSM) [p!o| , also known as the 
square-root measurement [|TÏ], 12]. Eldar and Forney [10] derive the optimal measurement for a 



pure-state ensemble in which the vectors \<j>ï) have a strong symmetry property called geomètric 



uniformity. In this case the vectors \<pi) are defined over a finite abelian group of unitary matrices 
and generated by a single generating vector; the optimal measurement again coincides with the 
LSM. Note, that a cyclic state set is a special case of a geometrically uniform state set. 

The LSM has many desirable properties [jn], 11, [Ï2|, [L|, ^, 14, and has therefore been 



proposed as a detection measurement in many settings (see e.g., |16|, |17], |18|). In Section y we 
derive a sufficient condition on the density operators for the LSM to minimize the probability of 
a detection error. For rank-one ensembles we show that the LSM minimizes the probability of a 
detection error if the probability of correctly detecting each of the states using the LSM is the 
same, regardless of the state transmitted. 

In Section |] we consider geometrically uniform (GU) state sets defined over a finite group of 



unitary matrices. In contrast to [10], the GU state sets we consider are not constrained to be 
rank-one state sets but rather can be mixed state sets, and the unitary group is not constrained to 
be abelian. We obtain a convenient characterization of the LSM and show that the LSM operators 
have the same symmetries as the original state set. We then show that for such GU state sets the 
probability of correctly detecting each of the states using the LSM is the same, so that for rank-one 
ensembles, the LSM minimizes the probability of a detection error. For an arbitrary GU ensemble, 
the optimal measurement operators are shown to be GU with the same generating group, and can 
be computed very efficiently in polynomial time. Furthermore, under a certain constraint on the 
generators, the LSM again minimizes the probability of a detection error. 

In Section || we consider the case in which the state set is generated by a group of unitary 
matrices using múltiple generators. Such a collection of states is referred to as a compound GU 
( CGU) state set |1S]. We obtain a convenient characterization of the LSM for CGU state sets, and 
show that the LSM vectors are themselves CGU. When the probability of correctly detecting each 
of the generators using the LSM is the same, we show that the probability of correctly detecting 
each of the states using the LSM is the same. Therefore, for rank-one CGU ensembles with this 
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property, the LSM minimizes the probability of a detection error. An interesting class of CGU 
state sets results when the set of generating vectors is itself GU, which we refer to as CGU state 
sets with GU generators. In the case in which the generating vectors are GU and generated by a 
group that commutes up to a phase factor with the CGU group, we show that the LSM vectors are 
also CGU with GU generators so that they are generated by a single generating vector. For such 
state sets, the probability of correctly detecting each of the states using the LSM is the same, so 
that for rank-one ensembles, the LSM minimizes the probability of a detection error. Finally we 
show that for arbitrary CGU state sets, the measurement operators minimizing the probability of 
a detection error are also CGU, and we propose an efhcient algorithm for computing the optimal 
generators. 

Before proceeding to the detailed development, in the next section we present our problem and 
summarize results from || pertaining to the conditions on the optimal measurement operators. 

2 Optimal Detection of Quantum States 

Assume that a quantum channel is prepared in a quantum state drawn from a collection of given 
states represented by density operators {pi, 1 < i < m} in an n-dimensional complex Hilbert space 
Tí. We assume without loss of generality that the eigenvectors of pi, 1 < i < m, collectively spanQ 
Tí so that m > n. Since pi is Hermitian and PSD, we can express pi as pi = (pi<j)* for some matrix 
(j>i, e.g., via the Cholesky or eigendecomposition of pi [Q. We refer to as a factor of p{. Note 
that the choice of <pi is not unique; if </>j is a factor of pi, then any matrix of the form = (ftiQi 
where Qi is an arbitrary matrix satisfying Q%Q* = I, is also a factor of pi. 

At the receiver, the constructed measurement comprises m PSD Hermitian measurement oper- 
ators {IL, 1 < i < m} on Tí that satisfy Y^iLi = ^> where / is the identity operator on Tí. We 

1 Otherwise we can transform the problem to a problem equivalent to the one considered in this paper by refor- 
mulating the problem on the subspace spanned by the eigenvectors of {pi, 1 < i < m}. 
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seek the measurement operators {Ilj, 1 < i < m} satisfying 



Ui > 0, 1 < i < m; 

m 

1=1 

that minimize the probability of a detection error, or equivalently, maximize the probability of 
correct detection. Given that the transmitted state is pj, the probability of correctly detecting 
the state using measurement operators {üj, 1 < i < m} is Tr(pjUj). Therefore, the probability of 
correct detection is given by 

m 

P d = Y,Pi^(Pi^i), (2) 
i=l 

where pi > is the prior probability of pi, with YliPi = 

It was shown in [||, ^] that a set of measurement operators {Ilj, 1 < i < m} minimizes the 
probability of a detection error for a state set {pi, 1 < i < m} with prior probabilities {pj, 1 < i < 
m} if and only if there exists an Hermitian X satisfying 

X > PiPi, 1 < i < m, (3) 

such that 

(X- Wí )ni = 0, l<i<m. (4) 
The matrix X can be determined as the solution to the problem 



min Tr(X) (5) 
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where B is the set of Hermitian operators on TC, subject to 



X > PíPí, 1 < i < m. (6) 

Except in some particular cases [[j], [7|, || ||, 10 1, obtaining a closed-form analytical expression 
for the optimal measurement operators directly from these necessary and sufficient conditions for 
optimality is a difficult and unsolved problem. Since (Jsj) is a (convex) semidefinite programming 
[pï| , 22, [2*|] problem, there are very emcient methods for solving (||). In particular, the optimal 



matrix X minimizing Tr(X) subject to @ can be computed in Matlab using the linear matrix 
inequality (LMI) Toolbox. A convenient interface for using the LMI toolbox is the Matlab package^] 
IQC/3 (see ]|] for further details). Once we determine X, the optimal measurement operators üj 
can be computed using (Q) and @) as described in ||]. 

A suboptimal measurement that has been employed as a detection measurement in many appli- 



cations and has many desirable properties is the LSM [10, 15]. Using the necessary and sufficient 
conditions (|ï|), (||) and (||), in Section || we derive a general condition under which the LSM is 
optimal, i.e., minimizes the probability of a detection error when distinguishing between possibly 
mixed quantum states. In Sections ||| and || we consider some special cases of mixed and pure 
state sets for which the LSM is optimal, and derive explicit formulas for the optimal measurement 
operators. 



2 This software was created by A. Megretski, C-Y. Kao, U. Jonsson and A. Rantzer and is available at 
http : //www .mit . edu/ cykao/home .html. 
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3 The LSM and the Optimal Measurement 

The LSM corresponding to a set of density operators {pi = (j)i4>*, 1 < i < m} with eigenvectors that 
collectively span TC and prior probabilities {pi, 1 < i < m} consists of the measurement operators 
{Sj = m/j*, 1 < i < m} where |§ [LO] 

m = (w*r 1/2 ^ = Tfr, (7) 



with 

r = (W)-V2. (8) 

Here \E' is the matrix of (block) columns ipi = y/pi4>i and (-) 1 / 2 is the unique Hermitian square root 
of the corresponding matrix. Note that since the eigenvectors of the {pi} collectively span 7í, the 
columns of the {tpi} also together span Tí, so ^>^i* is invertible. From ([?]), 

m / m \ 

Pipí = (^**)- 1 /2 ^ tilft (W)" 1 / 2 = (^^*)- 1 /2^^*(^^*)-l/2 = /; (g) 

i=l \i=l / 

so that the LSM operators defined by (0) satisfy (jï]). In the case in which the prior probabilities 
are all equal, 

m = ($$T 1/2 <fe (io) 

where $ is the matrix of (block) columns fa. 

Since the factors fa are not unique, the LSM factors pi are also not unique. In particular, if pi 
are the LSM factors corresponding to <pi, then the LSM factors corresponding to ^\ = faQi with 
QiQ* = I are p'i = PíQí- Therefore, although the LSM factors are not unique, the LSM operators 
Sj = pip* are unique. 

The LSM corresponding to a pure-state ensemble \(pi) consists of the measurement vectors 
\pi) = T\ipi), where = ^Jp~i\4>i)- It was shown in [|ï^] that for rank-one ensembles the LSM 
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vectors minimize the sum of the squared norms of the error vectors |ej) = — so that 
they are the measurement vectors that satisfy (||), and are closest in a squared error sense to the 
weighted state vectors = ^/pï\4>i)■ In the case in which the vectors \<f>i) are linearly independent 
so that n = m, (Q) implies that the vectors \fn) must be orthonormal, so that the LSM vectors are 
the closest orthonormal vectors to the vectors \tpi) in a least-squares sense, as illustrated in Fig. |]. 
A similar result was obtained for the LSM factor m corresponding to a mixed-state ensemble with 
factors fa Q. 




Figure 1: Example of the least-squares measurement (LSM). Since the vectors are linearly 
independent, the LSM vectors are orthonormal and minimize Y^i i e i\ e i) where \ei) = \ipi 



The LSM is equivalent to the square-root measurement M, 11, 16, 17, p^| , and has many 
desirable properties. Its construction is relatively simple; it can be determined directly from the 
given collection of states; it minimizes the probability of a detection error for pure-state ensembles 
that exhibit certain symmetries 0, [ÏÏJ ; it is "pretty good" when the states to be distinguished are 



equally likely and almost orthogonal [11|; it achieves a probability of error within a factor of two 



of the optimal probability of error [13]; and it is asymptotically optimal [12, 15]. Because of these 
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properties, the LSM has been proposed as a detection measurement in many applications (see e.g., 



m mm)- 

It turns out that in many cases of practical interest the LSM is optimal, i.e., it minimizes 
the probability of a detection error. From the necessary and sufficient conditions for optimality 
discussed in Section |2| it follows that the LSM minimizes the probability of a detection error if 
and only if the measurement operators LT, = pip* defined by ([?]) satisfy (Q) for some Hermitian X 
satisfying (^). A sufficient condition for optimality of the LSM is given in the following theorem, 
the proof of which is provided in the Appendix. 

Theorem 1. Let {pi = <pi<p*,l < i < m} denote a collection of quantum states with prior prob- 
abílítíes {pi,l < i < m}. Let {Sj = Pip\, 1 < i < m} with {pi = Tipi,l < i < m} denote 
the least-squares measurement (LSM) operators corresponding to {ipi = ^/pï4>i, 1 < i < m}, where 
T = (í'í'*) -1 / 2 and ^ is the matrix with block columns ipi. Then the LSM minimizes the probability 
of a detection error if for each i, p^ipi = ip^Tipi = al, where a is a constant independent of i. 

A similar result for the special case in which the density operators pi are rank-one operators of the 
form pi = \4>í){4>í\ and the vectors \(fii) are linearly independent was derived in |jÏ6| . 

Note that the condition ip*Tipi = al does not depend on the choice of factor <\>{. Indeed, if 
<fi'i = 4>iQi is another factor of p% with Qi satisfying QíQ* = I, and if *S>' is the matrix of block 
columns = y/pl^'i = y/Wi ( t ) iQii then it is easy to see that (tp^)*^'^'*)^ 1 ^ 2 ^ = al if and only if 
V>*TV>i = al. 

If the state p% = 4>i4>* is transmitted with prior probability pi, then the probability of correctly 
detecting the state using measurement operators £j = pifi* is p.;Tr(/z*</)jç!>* pi) = Tr(p*xpiip* pi). 
It follows that if the condition for optimality of Theorem |ï| is met, so that p*ipi = al, then the 
probability of correctly detecting each of the states pi using the LSM is a 2 , independent of i. 

For a pure-state ensemble consisting of states \4>í) with prior probabilities pi, the probability of 
correct detection of the ith state is given by \(pi\ipi}\ 2 . Since (pi\ipi} = (ipi\T\ipi) > for any set 



of weighted vectors \ipi), (pi\ipi) is constant for all i if and only if |(/Xj|^j)| 2 is constant for all i. 
Therefore, we may intèrpret the condition in Theorem || for pure-state ensembles as follows: The 
LSM is optimal for a set of states \4>i) with prior probabilities pi if the probability of detecting each 
one of the states using the LSM vectors is the same, regardless of the specific state chosen. 

In the remainder of the paper we use Theorem p] to derive the optimal measurement for mixed 
and pure state sets with certain symmetry properties. The symmetry properties we consider are 
quite general, and include many cases of practical interest. 

4 Geometrically Uniform State Sets 

In this section we consider the case in which the density operators pi are defined over a (not 
necessarily abelian) group of unitary matrices and are generated by a single generating matrix. Such 
a state set is cal·led geometrically uniform (GU) |24| ]. We first obtain a convenient characterization 
of the LSM in this case and then show that under a certain constraint on the generator, the LSM 
minimizes the probability of a detection error. In particular, for pure-state ensembles the LSM 
minimizes the probability of a detection error. 

Let Q = {Ui,l < i < m} be a finite group of m unitary matrices Ui. That is, Q contains the 
identity matrix /; if Q contains Ui, then it also contains its inverse U~ l = U*; and the product 
UiUj of any two elements of Q is in Q |p5| . 

A state set generated by Q using a single generating operator p is a set S = {pi = U- L pU* , Ui G G}- 
The group Q will be called the generating group of S. For concreteness we assume that U\ = I so 
that pi = p. Such a state set has strong symmetry properties and is called GU. For consistency 
with the symmetry of S, we will assume equiprobable prior probabilities on S. 

If the state set {pi,l < i < m} is GU, then we can always choose factors (pi of pi such that 
{4>i = Ui4>, Ui S G} where <j> is a factor of p, so that the factors (pi are also GU with generator (j). In 
the remainder of this section we explicitly assume that the factors are chosen to be GU. 
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We note that in |1C] a GU state set was defined for the case of rank-one ensembles. Furthermore, 
the generating group was assumed to be abelian. 

In the next section we derive the LSM operators for GU state sets and show that the LSM 
operators are also GU with the same generating group. We will see that this implies that when 
using the LSM, the probability of correct detection of each of the states in a GU state set is the 
same regardless of the particular state chosen. From Theorem |ï| it then follows that for pure-state 
ensembles, the LSM is optimal. 

4.1 The LSM for GU States 

To derive the LSM for a GU state set with generating group G, we first show that commutes 
with each of the matrices Ui G Q. Indeed, expressing <ï><ï>* as 

m m 



i=l i=l 



we have that for all j, 



i=i 

m 

= UjYWUiWUtUj 



i=l 

m 



üj E Uí<m>*Uí 



i=l 



= Uj<S><S>*, (12) 

since {U*Ui, 1 < i < m} is just a permutation of Q. 
If (JxJ)* commutes with Uj, then 

M = ($$*)-V2 ( 13 ) 
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also commutes with Uj for all j. Thus, from ([Ü]) the LSM operators are £j = \X{\x* with 



Mi = Mfa = MU t cJ) = UiMcj) = U il· L, (14) 

where 

H = Mcp = ($$*)- 1/2 0. (15) 

It follows that the LSM factors \ii are also GU with gener ating group Q and gener ator p, given by 
(|Ï5"p. Therefore, to compute the LSM factors for a GU state set all we need is to compute the 
generator p. The remaining measurement factors are then obtained by applying the group Q to p. 



A similar result was developed in [10] for rank-one ensembles in the case in which the group Q 



is abelian using the Fourier transform defined over Q. 



4.2 Optimality of the LSM 



We have seen that for a GU state set {pi = 4>i4>*, 1 < i < m} with equal prior probabilities 1/m 
and generating group Q = {Ui,l < i < m}, the LSM operators {Sj = /Xj/i*,l < i < m} are also 
GU with generating group Q. Therefore, 



A*i^i = -L·f J 'i ( l>i = -^=li*U*Uicj) = -^fi*(p, (16) 



where (f) and [i are the generators of the the state factors and the LSM factors, respectively. It 
follows that the probability of correct detection of each one of the states pi using the LSM is the 
same, regardless of the state transmitted. This then implies from Theorem [ï] that for pure-state 
GU ensembles the LSM is optimal. For a mixed-state ensemble, if the generator ep satisfies 

l· L* ( p = cP*(^*)- 1 / 2 cP = aI, (17) 
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for some a, then from Theorem [ï] the LSM minimizes the probability of a detection error. 

Note triat the condition p*(j> = al does not depend on the choice of generator ep. Indeed, if 
<$>' = 4>Q is another factor of p, then from ( |Ï5|) the generator of the LSM factors is p! = pQ so that 
fj,'*(f)' = al if and only if p*<p = al. 

4.3 Optimal Measurement for Arbitrary GU State Sets 

If the generator cf> does not satisfy (|ï~7|), then the LSM is no longer guaranteed to be optimal. 
Nonetheless, as we now show, the optimal measurement operators that minimize the probability of 
a detection error are GU with generating group Q. The corresponding generator can be computed 
very efficiently in polynomial time. 

Suppose that the optimal measurement operators that maximize 

m 

J({n í }) = ^Tr( /0i n i ), (18) 
i=i 

are IL, and let J = J({rL}) = £)I=i Tr(pjHj). Let r(j,i) be the mapping from X x X to X 
with X = {l,...,m}, defined by r(j,i) = k if Uf Ui = U^. Then the measurement operators 
LT^ = UjH r (j^U* for any 1 < j < m are also optimal. Indeed, since IL > and YsZLi n» = U 
U'i>0 and 

m / m \ 

E°í = íí i PK = ^ = 7 - (19) 

i=l \i=l / 

Finally, using the fact that pi = UipU* for some generator p, 

m mm 

«/({na) = Ys t *(p u ? u Au,í)U;u) = J2MpUk^kU k ) = ^Tr(par) = j. (20) 

1=1 fc=l i=l 

Since the measurement operators {IT^ = Ujllj.fj^Uf , 1 < i < m} are optimal for any j, it follows 
immediately that the measurement operators {II, = (1/m) Y^jL·i ^j^-r(j,i)Uf , 1 < i < ?n} are also 
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optimal. Indeed, it is immediate that Ilj satisfy ([[]). In addition, J({IIj}) = J({II^}) = J. Now, 

3=1 
^ m 

= -Y,UiUZn k u k u* 

k=l 

= UiÜU*, (21) 

where fi = (1/m) J2T=i U^ k U k . 

We therefore conclude that the optimal measurement operators can always be chosen to be GU 
with the same generating group Q as the original state set. Thus, to find the optimal measurement 
operators all we need is to find the optimal generator II. The remaining operators are obtained by 
applying the group Q to II. 

Since the optimal measurement operators satisfy Ilj = UiUU* and pi = UipU* , Tr(/9jIIj) = 
Tr(pII), so that the problem (pi) reduces to the maximization problem 



maxTr(plT), (22) 



where B is the set of Hermitian operators on 7i, subject to the constraints 



n > 0; 

m 

J^UilíU* = I. (23) 

i=l 

Since this problem is a (convex) semidefinite programming problem, the optimal II can be computed 
very efhciently in polynomial time within any desired accuracy pï| , p2] , |23|| , for example using the 
LMI toolbox on Matlab. Note that the problem of (^) and (|2^ ) has n 2 real unknowns and 2 
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constraints, in contrast with the original maximization problem (Q) and (^j which has mn 2 real 
unknowns and m + 1 constraints. 

We summarize our results regarding GU state sets in the following theorem: 

Theorem 2 (GU state sets). Let S = {pi = UipU*,Ui G G} be a geometrically uniform (GU) 
state set generated by a finite group Q of unitary matríces, where p = (fxj)* is an arbítrary generator, 
and let <í> be the matrix of columns fa = Ui<ft. Then the least-squares measurement (LSM) is given 
by the measurement operators Sj = pip\ with 

Pi = Uip 

where 

The LSM has the following properties: 

1. The measurement operators are GU with generating group Q; 

2. The probability of correctly detecting each of the states pi using the LSM is the same; 

3. If p*4> = (j)*{<S><5>*)- l / 2 (t) = al then the LSM minimizes the probability of a detection error; In 
particular, if <p = \<p) is a vector so that the state set is a rank-one ensemble, then the LSM 
minimizes the probability of a detection error. 

For an arbitrary generator <f> the optimal measurement operators that minimize the probability of 
a detection error are also GU with generating group Q and generator II that maximizes Tr{pTl) 
subject to n > and YT=i Ui^U* = I. 



15 



5 Compound Geometrically Uniform State Sets 

In this section, we consider state sets which consist of subsets that are GU, and are therefore 
referred to as compound geometrically uniform (CGU) |Ï9f| . As we show, the LSM operators are 
also CGU so that they can be computed using a set of generators. Under a certain condition on 
the generators, we also show that the optimal measurement associated with a CGU state set is 
equal to the LSM. For arbitrary CGU state sets, the optimal measurement is no longer equal to the 
LSM. Nonetheless, as we show, the optimal measurement operators are also CGU and we derive 
an efficient computational method for finding the optimal generators. 

A CGU state set is defined as a set of density operators S = {pn^ = 0ifc0* fe; 1 < i < 1,1 < k < r} 
such that pik = UiPkU*, where the matrices {Ui,l < i < 1} are unitary and form a group Ç, and 
the operators {pk, 1 < k < r} are the generators. For concreteness we assume that U\ = I so that 
Pik = Pk- We also assume equiprobable prior probabilities on S. 

If the state set {pik, 1 < i < 1,1 < k < r} is CGU, then we can always choose factors (f>ik of pu~ 
such that {4>ik = Ui4>k, 1 < i < 1} where (pk is a factor of pk, so that the factors cj)^ are also CGU 
with generators {4>k, 1 < k < r}. In the remainder of this section we explicitly assume that the 
factors are chosen to be CGU. 

A CGU state set is in general not GU. However, for every k, the matrices {4>ik, 1 5í i 5í an d 
the operators {pik, 1 < i < 1} are GU with generating group Q. 
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5.1 Example of a Compound Geometrically Uniform State Set 

An example of a CGU state set is illustrated in Fig. ||. In this example the state set is {pik = 
\<Pik)((f>ik\A <i,k<2} where {\4> ik ) = Ui\(f> k ),Ui € G}, Q = {h,U} with 




1 ^3 



(24) 



V3 -1 



and the generating vectors are 




1 




1 



(25) 



1 



-1 



The matrix U represents a reflection about the dashed line in Fig. |2[ Thus, the vector \<p2i) is 



obtained by reflecting the generator |</>ii) about this line, and similarly the vector \4>22) is obtained 
by reflecting the generator |çí>i2) about this line. 

As can be seen from the figure, the state set is not GU. In particular, there is no isometry that 
transforms |çí>ii) hito {(frw) while leaving the set invariant. However, the sets S\ = {\4>ii) , \(p2i)} 
and 52 = {|0i2)> 1022)} are both GU with generating group Q. 



We now derive the LSM for a CGU state set with equal prior probabilities. Let $ denote the matrix 
of columns (pik- Then for a CGU state set with generating group G, commutes with each of 
the matrices Ui £ Q. Indeed, expressing as 



5.2 The LSM for CGU State Sets 




(26) 
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I01l) 



I02l) 



1012) 



Figure 2: A compound geometrically uniform pure-state set. The state sets S\ = {|0n), |02i)} 
and £2 = {|0i2), 1022)} are both geometrically uniform (GU) with the same generating group; 
Both sets are invariant under a reflection about the dashed line. However, the combined set 
S = {|0n), |02i), |0i2), 1022)} is no longer GU. 



we have that for all j, 





= Uj<S><S>*, 



(27) 



since {U*Ui, 1 < i < 1} is just a permutation of Q. 

If <£<!>* commutes with Uj, then M = also commutes with Uj for all j. Thus, the 



LSM operators are Ejfc = 



with 



fi ik = M&k = MUi4>k = UiMcfik = UijJLk, 



(28) 
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where 



(29) 



Therefore the LSM factors are also CGU with generating group Q and generators fi^ given by 
(f29|). To compute the LSM factors all we need is to compute the generators fik- The remaining 
measurement factors are then obtained by applying the group Q to each of the generators. 
For the CGU state set of Fig. we have that 



1 
1 



(30) 



Therefore, the LSM vectors are {|//ifc) = Ui\fik),Ui £ Ç,l < i,k < 2} where Q = {Í2,Í7} with U 
given by (|24|), and from (|29| ) the generating vectors are 



A'l/ 



1 



M2) 



(31) 



The LSM vectors are depicted in Fig. ||. The vectors have the same symmetries as the state set 
of Fig. ^, with different generating vectors. Since in this example the generating vectors satisfy 
\Hk) = (1/V2)\4> k ), we have that = (l/y/2)\<f) ik ) for 1 < i,k < 2. 

5.3 CGU State Sets With GU Generators 

A special class of CGU state sets is CGU state sets with GU generators in which the generators 
{pk = 4 > k4 > *k-· 1 < k < r} and the factors (fij- are themselves GU. Specifically, {(ftk = Vk4>\ for some 
generator (fi, where the matrices {V^, 1 < A; < r} are unitary, and form a group Q. 

Suppose that U p and Vt commute up to a phase factor for all t and p so that U p Vt = VtU p e^^ p ' t ^ 
where 6(p,t) is an arbitrary phase function that may depend on the indices p and t. In this case 
we say that Q and Q commute up to a phase factor and that the corresponding state set is CGU 
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Figure 3: The least-squares measurement vectors associated with the compound geometrically 
uniform state set of Fig. ^. As can be seen from the figure, the measurement vectors have the same 
symmetries as the original state set. 



with commuting GU generators. (In the special case in which 6 = so that UiV k = V k Ui for all i, k, 
the resulting state set is GU Oli). Then for all p,t, 



$$*U p Vt = J2 Uí \J2 V ^* V k u*u p v t 
1=1 \k=l / 

= u p v t Y, v;u;uí (j2 v k <M>*vA u*u p v t 
ï=i Vfe=i / 

= u p v t v t *Ui [J2 Vk<M>*vn u*v t 
í=i Vfc=i / 

= u p v t Y, Ui (j2 v t *v k ^*v k *v)j u* 



i=l \k=l 

1=1 \fc=l / 

u p vm*. 



(32) 



The LSM factors fii k are then given by 



[i ik = M$ ik = MUiV k <P = UiV k M^ = UiV k fi, 



(33) 
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where fi = M(f). Thus even though the state set is not in general GU, the LSM factors can be 
computed using a single gener ator. 

Alternatively, we can express /ii k as hq. = Ui/j, k where the generators fi k are given by 



(34) 



From (34) it follows that the generators [i k are GU with generating group Q = {V k ,l < k < r} 
and generator p,. 

We conclude that for a CGU state set with commuting GU generators and generating group Q, 
the LSM vectors are also CGU with commuting GU generators and generating group Q. 



5.4 Example of a CGU State Set with Commuting GU Generators 

We now consider an example of a CGU state set with commuting GU generators. Consider the 
group G of l unitary matrices on C' where Ui = Z l ,l < i < l and Z is the matrix defined by 



Xl 




X2 


X2 




X3 


Xl-1 




Xl 


Xl 




Xl 



(35) 



Let Q be the group of r = l unitary matrices = B k , 1 < k < l where B is the diagonal matrix 
with diagonal elements e^ 2lTS ' 1 ,0 < s < l — 1. We can immediately verify that for this choice of 
Q and G, UiV^ = V^Uie^/ 1 . We therefore conclude that the LSM operators corresponding to the 
CGU state set S = {p ik = 4>ik<P* k } with {(j} ik = Ui4>k,Ui £ G} and {4> k = Vk<t>,Vk £ Q} for some 
generator cf>, are also CGU with commuting GU generators and can therefore be generated by a 
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single generator. 

As a special case, suppose that l = 2 so that Q consists of the matrices U\ = I2 and U2 = Z 
where 



Z 



1 

1 



(36) 



and Q consists of the matrices V\ = I2 and V2 = B where 



B 



1 

-1 



(37) 



Let the state set be S = {\4>ik} = UiVk\(j)), 1 < i,k < 2}, where 
normalized, 0\ + 0\ = 1. Then, 



[01 02]*- Since \4>) must be 



I0n> 



01 


, \hi) = 


02 




01 


, 1022) = 


-02 






, 1012) = 






_ _ 




. & . 




-02 _ 




01 _ 



(38) 



The LSM vectors are given by {\Hikj = UiV k \p), 1 < i,k < 2}, where = ($<&*) _1 / 2 |ç!>), and 



$ = 



A #2 01 -02 
02 01 -02 01 



(39) 



We can immediately verify that 





2{0\ + 01) 




2 










2(/? 2 + /3 2 2 ) 




2 



(40) 



Thus, the LSM vectors are 



(41) 
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In Fig. U we plot the state vectors given by (|3^) for the case in which \cj)) = (1/V5)[2 1]*. As can 
be seen from the figure, the state set is not GU. In particular, there is no isometry that transforms 
|0n) into |</>i2) while leaving the set invariant. Nonetheless, we have seen that the LSM vectors 
can be generated by a single generating vector |/i) = (l/y/2)\<p). 



1022) 


102!) 


s 

/ 

✓ 

✓ 




/ ✓ 

/ s 
1 s 
I s 
I s 
I ' 
/ * 

I * 


- |01l) 


S 

/ 

s 






s 

s 

s 

/ 

s 

s 




- 1012) 



Figure 4: A compound geometrically uniform state set with commuting geometrically uniform 
(GU) generators. The state sets S\ = {|0n), |02i)} arid £2 = {|0i2) 5 1022)} are both GU with the 
same generating group; Both sets are invariant under a reflection about the dashed line. The set of 
generators {|0n), |0i2)} is GU and is invariant under a reflection about the x-axis. The combined 
set S = {|0n), |02i)> 1012), 1022)} is no longer GU. Nonetheless, the LSM vectors are generated by 
a single generating vector and are given by \[Xik) = / V^)\<Pik) ■ 



Note, that in the example of Section 5A_ the CGU state set also has GU generators. Specifically, 
the set of generators {|0i), |02)} with |0i) = |0n) and |0a) = 1012) is invariant under a reflection 
about the x-axis: |02) = -B|0i) where B is given by (^7|). However, the group Q = {I2,B} of 
generators does not commute up to a phase with the generating group Q = {I2,U}, where U is 
given by ( p4| ) and represents a reflection about the dashed line in Fig. ^. This can be verified 
graphically from Fig. ^: Suppose we apply B to |0u) and then apply U. Then the resulting vector 
is equal to |022)- H on the other hand we first apply U to |0n) and then apply B, then the resulting 
vector is the reflection of |02i) about the x-axis, which is not related to |022) by a phase factor. 

Now, consider the state set in Fig. |j. In this case Q = {I2,B} and Q = {I2,Z} where B 
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represents a reflection about the x-axis and Z represents a reflection about the dashed line in 
Fig. |4|. We can immediately verify from the figure that applying Z and then B to any vector in the 
set results in a vector that is equal up to a minus sign to the vector that results from first applying 
B and then Z. For example, applying B to \4>u) and then applying Z results in \4>22)- If on the 
other hand we first apply Z to \4>u) and then apply B, then the resulting vector is the reflection 
of \4>2i) about the x-axis, which is equal to — | ^22) - 

5.5 Optimality of the LSM 

We have seen in the previous section that the LSM operators corresponding to a CGU state set with 
generating group Q = {Ui, 1 < i < /} is also CGU with the same generating group. In particular 
for each k, the sets S' k = {/J,ik, 1 < i < 1} and Sk = {4>ik, 1 < i < Z} are both GU with generating 
group Q. Therefore, 

fJ-ik&k = vlUiUi'l·k = V*k4>k, (42) 

which implies that the probability of correctly detecting each of the states in using the LSM is 
the same. It follows from Theorem || that if 

Vtfa = $%M<t) h = al, 1 < k < r, (43) 

then the LSM minimizes the probability of a detection error. 

Note that the condition fJ^(f>k = al does not depend on the choice of generator Indeed, if 
4>'k = 'PkQk is another factor of pk, then from (p9|) the generator of the LSM factors is /j,' k = pkQk 
so that p-'k^'k = a I if an d only if fit(f>k = 

In Section we showed that the LSM operators corresponding to a CGU state set with GU 
generators {4>k = Vk(f\ where 14 £ Q and Q and Q commute up to a phase factor, are also CGU 
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with GU generators generated by the same group Q and some generator /i. Therefore for all k, 

H%<i> k = fl*V k *V k 4> = ff</>, (44) 
so that the probability of correctly detecting each of the states <pi k is the same. If in addition, 

jfcj) = cj)*M(f) = al, (45) 
then combining (|42"|), (|44|) and ( |45"| ) with Theorem || we conclude that the LSM minimizes the 



probability of a detection error. In particular, for a rank-one ensemble, p,*(f> is a scalar so that (|45| ) 
is always satisfied. Therefore, for a rank-one CGU state set with commuting GU generators, the 
LSM minimizes the probability of a detection error. 

5.6 Optimal Measurement for Arbitrary CGU State Sets 

If the generators <j) k do not satisfy (|43|), then the LSM is no longer guaranteed to be optimal. 
Nonetheless, as we now show, the optimal measurement operators that minimize the probability of 
a detection error are CGU with generating group Q. The corresponding generators can be computed 
very efficiently in polynomial time within any desired accuracy. 

Suppose that the optimal measurement operators that maximize 

l r 

J({U ik }) =J2J2 MPi^ik), (46) 

i=l k=l 

are Ü ik , and let J = J({Ü ik }) = Ya=í 121=1 ^{Pik^hk)- Let r(j,i) be the mapping from Ixi 
to I with X = {1, ... ,1}, defined by r(j,i) = s if UJUi = U s . Then the measurement operators 
H' ik = UjH r (j ti -) k U* for any 1 < j < l are also optimal. Indeed, since Hi k > and Ylik^ik = I> 
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W ik > and 



l r ( l r \ 

i=i fc=i \ï=i fc=i / 



(47) 

Finally, using the fact that pi k = UíPkU* for some generators p k , 

l r l r l r 

J({%}) = Y J Y, Tl (P* U ? U firU^kU*U l ) = ^2^2Tr( Pk U;U sk U s ) = ^^Tr( m íl ife ) = J.(48) 

i=l fc=l S=l fc=l 1=1 fc=l 

Since the measurement operators {ïl' ik = UjH r ^j^ k U*,l < i < 1,1 < k < r} are optimal for 
any j, it follows immediately that the measurement operators {Hi k = (1/0 Ej=i ^j^-r(j,i)kUj , 1 < 
i < 1 < fe < ^} are also optimal. Indeed, it is immediate that Iljfc satisfy (]!]). In addition, 
J({ÏÏ ïfc }) = ,!({%}) = J. Now, 



n^fc = j ^ Uj^-r(j,i)kUj 

3=1 

1 í 

s=ï 



= UiU k U*, (49) 

where fi fc = (1/0 EÍ=l U;n a kU 8 . 

We therefore conclude that the optimal measurement operators can always be chosen to be CGU 
with the same generating group Q as the original state set. Thus, to find the optimal measurement 
operators all we need is to find the optimal generators {H k , 1 < k < r} . The remaining operators 
are obtained by applying the group Q to each of the generators. 

Since the optimal measurement operators satisfy Iljfc = UiU k U* and pi k = Uip k U* , Tr(pi k Hi k ) = 
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Tr (piJIk), so that the problem (|2|) reduces to the maximization problem 

r 

max^Tr(p fc II fc ), (50) 



lfct - k=l 



subject to the constraints 



n fc > 0, 1 < k < r 



l r 



ï=l fc=l 



Since this problem is a (convex) semidefinite programming problem, the optimal generators IT^ 
can be computed very efhciently in polynomial time within any desired accuracy [21, 22, 23], for 



example using the LMI toolbox on Matlab. Note that the problem of (50) and (51) has rn 2 real 
unknowns and r + 1 constraints, in contrast with the original maximization problem (||) and (|^) 
which has Irn 2 real unknowns and Ir + 1 constraints. 

We summarize our results regarding CGU state sets in the following theorem: 

Theorem 3 (CGU state sets). Let S = {pit = UiPkU* ,Ui £ G, 1 < k < r} be a compound geo- 
metrically uniform (CGU) state set generated by a finite group Q of unitary matrices and generators 
{Pk = ^k^ti 1 ^ k < r}, and let <3> be the matrix of columns <^ = Ui4>h- Then the least-squares 
measurement (LSM) is given by the measurement operators Ej = fiip* with 

Pik = UiPk 

where 

p k = ($$T 1/2 <^- 

The LSM has the following properties: 
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1. The measurement operators are CGU with generatíng group Q; 

2. The probabílity of correctly detecting each of the states for fixed k using the LSM is the 
same; 

3. If /i*0 fc = 0*($$*)-V20 fc = al forl <k <r then the LSM minimizes the probability of a 
detection error. 

If in addition the generators {(fik = Vk<fi-, l < k < r} are geometrically uniform with UiV^ = 
VkUieMW for all i, k, then 

1. = UiVkp- where fi = (<J><ï>*) -1 / 2 çí> so that the LSM operators are CGU with geometrically 
uniform generators; 

2. The probability of correctly detecting each of the states 4>ik using the LSM is the same; 

3. Ifp,*(f) = 0*($$*)-V20 = al then the LSM minimizes the probability of a detection error. In 
particular, if (fi = \cfi) is a vector so that the state set is a rank-one ensemble, then the LSM 
minimizes the probability of a detection error. 

For arbitrary CGU state sets the optimal measurement operators that minimize the probability of a 
detection error are CGU with generatíng group Q and generators üj, that maximize Ylk=i Tr(pk^ík) 
subject to Hk > 0, 1 < k < r and ^ k UiHkU* = I. 

6 Conclusion 

In this paper we considered the optimal measurement operators that minimize the probability 
of a detection error when distinguishing between a collection of mixed quantum states. We first 
derived a general condition under which the LSM minimizes the probability of a detection error. 
We then considered state sets with a broad class of symmetry properties for which the LSM is 
optimal. Specifically, we showed that for GU state sets and for CGU state sets with generators 
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that satisfy certain constraints, the LSM is optimal. We also showed that for arbitrary GU and 
CGU state sets, the optimal measurement operators have the same symmetries as the original 
state sets. Therefore, to compute the optimal measurement operators, we need only to compute 
the corresponding generators. As we showed, the generators can be computed very eniciently in 
polynomial time within any desired accuracy by solving a semidefinite programming problem. 



Appendix 

Proof of Theorem [j] 

In this appendix we prové Theorem [ï} Specifically, we show that for a set of states pi = (pifi* with 
prior probabilities pi, if fi*ipi = al, where m = (í'í'*) -1 / 2 ?/^ are the LSM factors and ifii = ■ s /pl4>i, 
then there exists an Hermitian X such that 



1 < i < m: 



{x - ^Dnm* 



0, 1 < i < m. 



(52) 



Let X be the symmetric matrix defined by X = aW 1 ^ 2 where W = Since al 

^W- l / 2 ipj = tfW-WW-Vty, it follows that 



al > W-^ipjrfW- 1 /*. 



(53) 



Multiplying both sides of by W 1 ^ we have 



aW 1/2 > fytà, 



(54) 



which verifies that the conditions (52) are satisfied. 
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Next, 

(X - Mïïm = a^*) 1 ' 2 ^*)- 1 ' 2 ^ - a^i = 0, (55) 



so that (p2|) is also satisfied. 
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